Dynamical Mean Field Study of The Dirac Liquid 
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Renormalization is one of the basic notions of condensed matter physics. Based on the concept of renor- 
mahzation, the Landau's Fermi liquid theory has been able to explain, why despite the presence of Coulomb 
interactions, the free electron theory works so well for simple metals with extended Fermi surface (FS). The 
recent synthesis of graphene has provided the condensed matter physicists with a low energy laboratory of 
Dirac fermions where instead of a FS, one has two Fermi points. Many exciting phenomena in graphene can be 
successfully interpreted in terms of free Dirac electrons. In this paper, employing dynamical mean field theory 
(DMFT), we show that an interacting Dirac sea is essentially an effective free Dirac theory. This observation 
suggests the notion of Dirac liquid as a fixed point of interacting 2+1 dimensional Dirac fermions. We find one 
more fixed point at strong interactions describing a Mott insulating state, and address the nature of semi-metal 
to insulator (SMIT) transition in this system. 

PACS numbers: 71.1().Fd, 73.43.Nq 



INTRODUCTION 

Dirac theory of electrons was formulated in 1928 to de- 
scribe the relativistic motion of electron waves This the- 
ory not only is consistent with the spin h/2 of electron, but 
also obeys the correct relativistic covariance. Now days, ac- 
celerators and neutron stars are not the only places to search 
for Dirac fermions. Advances in science and technology has 
enabled physicists to realize Dirac fermions in energy scales 
as low as leV, or even lower in solid state physics. The nodal 
qasiparticles of d-wave cuprate superconductors being one ex- 
ample in superconducting state of matter [2;]. 

Recently^raphene, a single atomic layer of graphite, was 
fabricated |3|], which has initiated a rapidly growing research 
activity in condensed matter physics |4J. The low-energy elec- 
tronic structure of this system is approximately described by a 
2+1 dimensional Dirac theory IH, which enjoys a chiral sym- 
metry islSl- Therefore graphene, the unrolled carbon nano- 
tube; provides the condensed matter with a laboratory for the 
relativistic fermions on the table top [Tt] . 

There is still another road to search for Dirac fermions in 
material physics laboratory. Recent advances in ultra cooling 
and atom trap methods |Ql have elevated this technology to a 
chip based level Jgt]. This provides a unique opportunity for 
tuning the interaction parameters in microscopic models em- 
ploying the so called Feshbach resonance llioll . These moti- 
vate the study of strongly interacting fermions on honeycomb 
lattice in parameters regimes, much beyond what can be cur- 
rently realized in graphene, or high Tc cuprates. 

The three-fold coordination for hopping of fermions on 
a honeycomb lattice is responsible for the relativistic low- 
energy theory which strictly speaking describes a semi- 
metallic state of a Dirac sea; rather than a metallic state of a 
Fermi sea. This semi-metal is described by a pseudo-gap in 
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FIG. 1: Lattice structure of a honeycomb lattice as a superposi- 
tion of two simple monoclinic Bravais lattices. The C-C distance 
is a ~ L42A, such that the lattice constant becomes \/^a and the 
basis vectors are given by ai = 3a/2ex + V3a/2ey, a.2 = y/Sae^. 



the density of states (DOS) shown in Fig. |2] Such feature in 
DOS is also relevant to the nodal quasi particles of high Tc 
cuprate superconductors . 

While the single particle Dirac theory seems to work well in 
graphene, it is important to search for interesting many-body 
effects on honeycomb lattice. For example the possi bility of 
/-wave or d + id- wave siiperconductivity III , 12, 13, 3, as 
well as CDW instability Ml ill on honeycomb lattice has been 
discussed. Perturbative renormalization group studies 1 15 . 3 
has indicated that the long range part of the Coulomb interac- 
tion is irrelevant, and gives rises to a non interacting fixed 
point. In this paper we focus on the short range part of the 
interaction, employing a non-perturbative method of DMFT. 
This theory has been very successful in addressing the ques- 
tion of metal to insulator transitions UtIi . According to the 
prediction of DMFT, the onset of transition to Mott insulating 
phase is accompanied by formation of the so called Kondo res- 
onance, the spectral weight of which characterizes the quasi- 
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FIG. 2: Density of states for tight binding fermions on honeycomb 
{D = 2) and hyper honeycomb {D = oo) lattices. In the case of 
undoped graphene, due to spin degeneracy the valence band (mirror 
image of this figure) is completely filled. For spin-less fermions on 
optical lattices at quarter filling also the chemical potential is at a; = 
0. Energies are in units of hopping amplitude t. 



particle weight Z of the underlying Fermi sea. In the case of 
a half-filled Dirac sea, instead of an extended Fermi surface, 
one has to deal with two Fermi points at the so called K points 
of the Brillouin zone. Due to the cone like dispersion of the 
conduction and valence band near these points, there will be 
no quasi particle state at the Fermi level 14]]. Hence the quasi- 
particle weight Z of Fermi liquids can not be used to discuss 
the transition to Mott insulating regime. 

In this work we employ the DMFT approximation to study 
the SMIT of Dirac liquids in two dimensional honeycomb lat- 
tice. The picture which emerges from this study is that; for 
weak to moderate interaction strengths, the Dirac sea state re- 
mains stable against local many body interactions. The sole 
role of interactions would be to renormalize the Dirac quasi 
particles. Such robustness of the Dirac fermions against many 
body interactions has been observed in various measurements 
in graphene jlSll . For strong enough interactions, a Mott in- 
sulating state is stabilized. We study some simpler model 
density of states which mimic the true 2D DOS of graphene 
Eq. (|3]l. We find the value of Uc ^ 13. 3i, which is not sensi- 
tive to the details of DOS. 



MODEL AND METHOD 



We take the tight-binding electrons on honeycomb lat- 
tice 101, which give rise to Dirac spectrum near the K points 
of the Brillouin Zone (BZ), and add a short ranged Coulomb 
interaction of Hubbard type to it: 



h.c. 



(1) 

where the hoping amplitude for electrons in graphene is 
t ~ 2.7eV, while for the atoms trapped in optical lattices it 



is typically on the scale of t ~ pK ~ 10~^'^eV. Through 
the paper when the units are not specified, the units in which 
t = h = 1 is implied. Also we take the atomic separation to 
be 1. This model is manifestly particIe-hoIe symmetric and 
the Hartree term at half filling {{rij^) = 1/2) vanishes. 
In the absence of interaction, the Hamiltonian becomes 



Ho 







-/(k) 
-/*(k) 



(2) 



_ „ik.di 



ik.ds 



where k = {k^,ky), and /(k) = 
(Fig. [U. The eigenvalues of this matrix are Sk = 
When the above complex function /(k) linearized around the 
Fermi points K, K' of the BZ, gives rise to the Hamiltonian 
H = zLvpa.'k. The DOS for this non interacting system is 
given by lll9ll . 




where 
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Here _FY7r/2,x) is the complete elliptic integral of first 
kind 12011 . At low energies where the dispersion becomes 
ek = ivpk, we have linear energy dependence in the pseudo- 
gap shaped DOS (Fig.Il), 



pis) = 2Trvf \e\ 



(6) 



where hvp = 3ta/2 is the bare Fermi velocity of quasiparti- 
cles at the Fermi points. In Fig.|2]we have plotted the DOS for 
two and infinite dimensional (hyper)honeycomb lattice ll2lll . 

In DMFT, one starts with a local free propagator gQa{uj), 
where the Greek indices a, [3, etc. correspond to sub-lattices 
A,B. Next and most important part of the DMFT consists in 
employing an impurity solver to obtain the local self-energy 
as a functional of g^a- Note that, within the single- 
site DMFT approximation the self-energy is purely local llTIl 
which means the off-diagonal self-energies vanish and it can 
be described by the diagonal matrix elements Yjaioj). We use 
the iterated perturbation theory to solve the impurity problem: 



(7) 



Now using this self-energy we construct the Green's function 
of the honeycomb lattice which can be written as 



_(QaO^,uj) -/(k) 
I -/*(k) CB(k,c.) 



(8) 



where k belong to the first BZ of of sub-lattice, and 
CQ(k, Lj) — uj — £y + p — Tja{Lo)- The half filled band which 



FIG. 3: (Color online) Changes in the band structure of Dirac 
Fermions on honeycomb lattice as a function of Hubbard U. For 
values of U larger than the bandwidth W = 6t, the upper and lower 
(not shown in this figure) Hubbard bands set in. 



gives rise to a Dirac sea, due to the bipartite nature of the hon- 
eycomb lattice corresponds to ^ = 0. Such half-filling con- 
dition corresponds to undoped graphene, whereas in optical 
lattices with one specie atom, it corresponds to half atom per 
site. Note that the particle-hole symmetry can also be used 
which implies = — ti)). To connect the lattice 

Green's function to local one, we need to project dH) onto sub- 
lattice A{B), which reads 

9aiu^) - E . / 2 = / ds-^^ (9) 

where p(e) is the DOS of massless Dirac fermions, Eq. (O. To 
close the set of equations, we only need to append the Dyson 
equation 

.9a'H=5o"aH-SaH (10) 

We solve the set of equations (O, ^ and ( fTOl i self- 
consistently. 

RESULTS 

The Hilbert transformation in equation (|9]l is the only place 
where the band structure (DOS) of system enters the DMFT 
machinery. Using the DOS of realistic graphene given in 
Eq. (|3]l we obtain the interacting DOS for various values of 
[/ < J7c ~ 13.3t shown in Fig. [3] As can be seen, upon 
increasing U, the spectral weight is transferred to higher en- 
ergies, accompanied by a increase in the slope of the DOS. 
Interpreting the results in terms of Eq. (|6]l one sees that (i) the 
Dirac nature of the spectrum is preserved, manifested in the 
linear DOS and the pseudo-gap structure around the charge 
neutrality point is preserved for all U < Uc, but with a re- 
duced Fermi velocity vp replacing the non-interacting Fermi 



FIG. 4: (Color online) DOS versus energy. Energy is in units of hop- 
ping parameter t. (a), (c): Renormalization of Dirac sea by Hubbard 
interaction U. (b), (d): SMIT for Dirac fermions. Upper two panels 
(a,b) correspond to h — 0, while the lower two panels (c,d) corre- 
spond to = 1 in the toy DOS introduced in the text. The values of 
U are in units of t. 



velocity vp = 3/2. (ii) As there are no quasiparticle states 
at the Fermi level for the Dirac fermions, the SMIT is not ac- 
companied by a Kondo resonance at the Fermi surface, (iii) 
The (logarithmic) van Hove singularity of the band structure 
corresponding to saddle point at the M point of the Brillouin 
zone persists. 

For the DOS (O the Hilbert transformation must be done 
with a numerical quadrature with typically Na = 2000 ab- 
scissas for a Simpson quadrature to achieve the double pre- 
cision accuracy needed in typical self-consistency problems. 
We simplify the DOS of graphene with a parameter dependent 
DOS which mimics the essential features of the graphene band 
structure, but at the same time allows for analytic evaluation 
of the Hilbert transform (|9]l which saves a lot of time: 

r 2nvp^\e\ \e\ < 1 
Phis) = hdi\e\ - 1) + i 27rvf K |£| < 3 (11) 
[ otherwise 

This model DOS captures the linear dispersion around the 
charge neutrality point and a sharp singularity of strength h at 
|£| = ±1. Solving the DMFT equations with DOS (dB in the 
projection formula (|9]l we obtain figure ID Panels (a), (b) cor- 
respond to /i = and (c), (d) denote h ~ 1. We see that this 
model DOS produces the essential features (i)-(iii) above. For 
small to intermediate values of U one can see the renormal- 
ization of the slope in (a), (c). Comparison of (a), (c) shows 
that this feature is qualitatively independent of the presence or 
lack of a singularity. Singularity places initially at |e| = ±1 
moves to lower energies as one increases U. Panels (b), (d) 
show the SMIT for h = 0, h = 1, respectively. For h = 0, 
we obtain J/^ ~ 13. 3t, while for h = 1, the critical value is 
slightly lower. Right at the critical point, due to critical slow- 
ing down, it is extremely difficult to obtain convergence. The 
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FIG. 5: (Color online) Renormalization of the Dirac fermions veloc- 
ity «F by the Hubbard U term. 



reason is that for U = , the DOS slope diverges (meaning 
— * 0), and the pseudo-gap tends to closes. 

In Fig. |5]we show the renormalized Fermi velocity vp for 
the model DOS corresponding to /t, = as a function of U. 
The Fermi vp of Dirac fermions interpolates to zero at the 
quantum critical point (QCP) Uc ~ 13. 3t. Unlike the MIT 
in metallic systems with extended Fermi surface, there is no 
Kondo resonance corresponding to quasiparticle at the Fermi 
level. Therefore SMIT can not be described in terms of the 
spectral weight Z of such resonant state. Instead the vp in 
Fig. Scan be identified as an order parameter characterizing 
a Dirac liquid state. Beyond a QCP at Uc ~ 13. 3t, a Mott 
insulating state appears. 

The semi-metal to insulator transition in Z) = oo hon- 
eycomb lattice (Dotted line in Fig. |2]i has been previously 
studied by Santoro and coworkers [21], where they focused 
on the zero and finite temperature magnetic phase transi- 
tions. They find that at zero temperature, at a critical value of 
« 2.3W there is a phase transition from a paramagnetic 
semi-metal to an anti-ferromagnetic Mott insulator, where W 
is the band width. If one takes this as a ad hoc critical value for 
D = 2 honeycomb lattice with W = 6t, one expects a criti- 
cal value of w 13. 8i in 2D honeycomb lattice, which is 
very close to the value 13. 3< we find here. Sorella and Tosatti 
approached the same problem by quantum Monte Carlo meth- 
ods 12211 . They found a zero temperature SMIT between a non 
magnetic semi-metal and an anti-ferromagnetic insulator at a 
critical value of [/Q^c ^ 4,^t. The Brinkman-Rice analysis 
of the SMIT within the Gutzwiller approximation gives a crit- 
ical value of U^^ « 12. 8t 1231], which is again rather close 
to the value we obtain. The corresponding mean field value is 
~ 2.23t. Such a significant difference between the mean 
field and more accurate methods indicates the importance of 
quantum fluctuations in SMIT. 



SUMMARY AND DISCUSSION 

In this paper we studied the zero temperature phase tran- 
sition of Dirac fermions on a honeycomb lattice. Our results 
can be applied to atoms in honeycomb optical lattices, as well 
as electrons in undoped graphene. We studied the nature of 
SMIT in this system. We found no quasiparticle resonance 
state. The renormalized Fermi velocity vp was suggested as 
an order parameter for a Dirac liquid separated from a Mott 
insulator by a QCP at J7c ~ 13. 3t. The recent renormalization 



in Vp seen in ARPES measurements 11 1811 can not be solely 
described in terms of electron-phonon or electron-plasmon in- 
teractions of a doped graphene. Such renormalization aris- 
ing from Hubbard term may also contribute to those seen in 
ARPES measurements. 

Others have also reported on a quantum critical point be- 
havior of Dirac fermions using large N expansion, beyond 
which an insulating phase emerges ll24l] . The renormaliza- 



tion group flow equations for quartic perturbations ||25|] indi- 
cate possible insulating phase at strong couplings. Foster and 
Aleiner find logarithmic divergences using large N renormal- 
ization group approach iflill . They found that the long range 
part of the Coulomb interaction enhances the short range part 
of the interaction. Our method addresses the role of short 
range part in a non-perturbative DMFT sense, which therefore 
can be regarded as a complementary analysis to their's. 

Gonzalez and coworkers find within renormalization 
scheme that the Coulomb interactions drive the system to a 
non-interacting fixed point fJ3\- Their finding agrees with the 
Dirac Liquid fixed point in our analysis. However our non- 
perturbative DMFT analysis indicates another Mott insulat- 
ing fixed point in addition to the Dirac liquid fixed point of 
Ref. 11511 in agreement with our exact and non-perturbative 
analysis. 

In graphene samples J7 « 6 eV, which is not strongly corre- 
lated, places them at U/t ^ 2.2, far from the Mott insulating 
phase. However, one can describe the graphene starting from a 



Mott insulating resonance valence bond (RVB) l26l 12711 point. 
In such a description one needs to allow for charge fluctua- 
tions on top of an RVB ground state 11311 . This provides an 
opportunity for h igh temperature superconductivity in ^ 20% 
doped graphene iflil. [l4ll . 
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